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We study the intermittent dynamics and the fluctuations of the dynamic correlation function of 
a simple aging system. Given its size L and its coherence length ^, the system can be divided into 
A' independent subsystems, where N — (j)''', and d is the dimension of space. Each of them is 
considered as an aging subsystem which evolves according to an activated dynamics between energy 
levels. We compute analytically the distribution of trapping times for the global system, which can 
take power-law, stretched-exponential or exponential forms according to the values of A'' and the 
regime of times considered. An effective number of subsystems at age t^, Neff{tw), can be defined, 
which decreases as increases, as well as an effective coherence length, ^(tii,) ~ t™"^'^'', where 
H < 1 characterizes the trapping times distribution of a single subsystem. We also compute the 
probability distribution functions of the time intervals between large decorrelations, which exhibit 
different power-law behaviours as increases (or A' decreases), and which should be accessible 
experimentally. Finally, we calculate the probability distribution function of the two-time correlator. 
We show that in a phenomenological approach, where A'^ is replaced by the effective number of 
subsystems N^ffit^), the same qualitative behaviour as in experiments and simulations of several 
glassy systems can be obtained. 
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I. INTRODUCTION 

The dynamics of glassy materials such as spin-glasses, 
structural glasses or amorphous soft materials like gels, 
pastes or foams has been a subject of considerable study 

E II 

Considerable effort has been made in order 
to understand and quantify the out-of-equilibrium char- 
acter of their temporal relaxation, in particular the ab- 
sence of time translational invariance (aging), through 
the study of dynamic correlation functions. On the the- 
oretical side, global dynamic correlations have been de- 
scribed at a mean-field level in disordered systems such 
as spin-glasses (Jj , or at a phenomenological level in mod- 
els such as the Random Energy Model 0, Q . Both ap- 
' proaches neglect all spatial properties of the system, and 
are therefore likely to miss any spatial correlations that 
\ arise during the dynamics. Alternatively, a phenomeno- 
logical picture, the droplet model, has been proposed, 
that focuses on the spatial properties as a key to under- 
stand the slow dynamics and the critical properties of 
\ spin-glasses 

In recent years, interest in the spatial properties of 
\ glassy systems has been growing. The size of excita- 
I tions in finite-dimensional spin-glasses has been stud- 
ied numerically |^. In simulations of kinetically con- 
strained glassy systems [lOl (llj , and of supercooled liq- 
uids cooperativity lengths have been identified, 
and related to the presence of heterogeneities in the dy- 
namics 

More generally, a lot of physical questions remain; a 
crucial one being: in what manner, at a microscopic level, 
does a glassy system evolve, both spatially and tempo- 
rally? What are the spatial configurations of the typical 



rearrangements experienced by a glassy system during 
its relaxation and how are they affected by aging? Are 
there any common relaxation mechanisms of glassy sys- 
tems, though there also exist specificities to given mate- 
rials? Recently, new results have been obtained in this 
direction. Focusing on the spatial aspect of glassy relax- 
ation, cooperative rearrangements events have been ev- 
idenced, both experimentally 0, 0| and through com- 
puter simulations [l^ . stimulating new research on the 
challenging question of coherence lengths and coopera- 
tivity in glasses. On the temporal side, beautiful experi- 
ments have shown evidence of temporal intermittency in 
colloidal gels and micellar polycrystals 0|, and in poly- 
carbonate glasses . It seems now well established that 
in glasses and gels, relaxation takes place in a discontin- 
uous way, involving sudden rearrangements followed by 
periods of arrest where almost nothing happens. The 
precise experimental determination of the distribution of 
time lags between rearranging events will give insight into 
the characterictic "trapping times" of the system. In ex- 
periments on glasses, this distribution seems to be close 
to a power-law -which is consistent with a trap model 
with an exponential distribution of energies j^-, whereas 
it has been found in simulations of supercooled glasses to 
correspond to a model of traps with either an exponen- 
tiel distribution of energies, or a gaussian distribution of 
energies 0| ^2 ■ Therefore, quantities of interest are not 
just average quantities, but also fluctuations, and in par- 
ticular the full probability distributions of correlations. 
The study of fluctuations in glassy systems may contain 
subtle information, as was already realized by Israeloff 
and Weissman who analyzed carefully mesoscopic 
noise in spin-glasses in an attempt to discriminate be- 
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tween a model of droplets, and a scenario of hierarchical 
dynamics. 

Recently, probability distributions of two-time corre- 
lation functions in gels and glasses have been shown to 
exhibit a non-Gaussian behaviour [l9l Eo|. These non- 
Gaussian features have been found also in numerical sim- 
ulations of disordered systems and kinetically constrained 
models (2^|25|,|2^] . They have been tentatively explained 
using the analogy between glassy dynamics and critical 
dynamics , for which universal, non-Gaussian features 
can be expected |2^ . 

In this study, we will not advocate any similarity with 
critical dynamics, but we will rather try to deduce the 
non-Gaussian behaviour merely from finite-size effects in 
a simple out-of-equilibrium glassy model. 

We will consider a system that can be divided into N 
independent subsystems. Each subsystem is supposed 
to represent an independent model of glassy relaxation 
between energy traps. Such a model of traps has been 
studied extensively . We will show that the super- 

position of the N subsystems will have the same aver- 
age dynamical correlation as one individual subsystem; 
however, its probability distribution function will depend 
strongly on A'^; in the limit N oo, one has to recover 
the Gaussian distribution, according to the central limit 
theorem. Moreover, the distribution of time intervals be- 
tween relaxation events or between decorrelations will 
depend crucially on the number of subsystems, i.e on the 
value of the internal coherence length. 

The paper is organized as follows. In Section II, we 
introduce the model and recall the main results of ^ . In 
Section III, we calculate the distribution of time inter- 
vals between successive events in the whole system. In 
Section IV, we calculate the distribution of time intervals 
between successive decorrelations in the system, the same 
way they can be measured in Time Resolved Correlation 
experiments on soft glassy materials Section V is 

devoted to the probability distribution of the two-time 
correlation function. Finally, in Section VI, we summa- 
rize and discuss our results. 



II. DEFINITION OF THE MODEL 

Let us consider a simple model of dynamics of a system 
between energy levels. In a generic disordered system, it 
is reasonable to assume that low-lying energy levels are 
exponentially distributed: this is the case for example 
for the lowest energy levels (non-extensive corrections to 
the ground state energy) of the Random Energy model 
(REM) Uli^ or in spin-glasses Experimental de- 
termination of energy barriers in low-temperature glasses 
also seem to support this exponential distribution |3 j |. 

More precisely, let us call E an energy barrier, which 
is the difference between a reference energy level which 
is taken as the origin for the energies, and a negative 



energy level. E is hence a positive quantity. We choose 
the distribution of barriers as p{E) = -g^e'^^^"- 

Changes in configurations in a disordered or glassy sys- 
tem are often attributed to thermally activated events 
over energy barriers [21I , although other mecanisms 
exist, such as kinetically constrained models which 
do not require an energy landscape, and are able to re- 
produce some of the features of glassy materials. In the 
language of activated events, a trapping time r corre- 
sponding to a barrier E can be defined as r = Toe^/'^^"^, 
where tq is a microscopic time scale, ks the Boltzmann 
constant, and T the temperature. For an exponential dis- 
tribution of barriers, the distribution of trapping times 

is equal to tp{T) — fJ-r^ii^, where ^ = ^g^. 

The dynamics of a system evolving in such an energy 
landscape has been studied by different authors, accord- 
ing to which choice of transition rates between energy 
levels is made [^[i^ls^li^l- In the following, we will con- 
sider only the case where the transition rate from barrier 
E to any other barrierE" is W{E E') = ^e'^/^^'^, 
1^ , which means that the escape from the initial trap 
is the limiting process, whatever the destination. Note 
that this family of models does not include any kind of 
spatial structure, since energy is not related here to spa- 
tial configurations, and that it is mean-field in nature 
since transitions to all levels are allowed with the same 
probability. Extensions of these models to finite dimen- 
sionalities have however been attempted in |3 US ■ We 
will see in the following that the superposition of several 
of such systems can actually introduce (though rather 
artificially) a relevant lengthscale. 

Dynamicalproperties of the model have been studied 
in detail in [6j. In particular, when /i < 1, the model 
exhibits aging (absence of time translation invariance of 
the correlations), whereas it is time-translationally in- 
variant for /i > 1. In this paper, we will focus on the 
case /i < 1. For a given trajectory of the system, the 
two-time correlation function C{tw , tui + 7") is defined as: 
C (iu) , tio + t) = 1 if the system has remained in the same 
energy trap between and tyj + r, C(tu,, -f r) = if 
between and -\- t, the system has left the energy 
trap it was in at i^. 

Averaging over all barrier configurations (which we de- 
note by < .. >), one obtains the average two-time corre- 
lation function: 

n(tu.,i«, +t) =< C(tu,,i„ + r) > 

It was shown in that for large t^,, this function is 
given by the following formula: 



Ii-{tw,tyj + t) 



sm TT/i 



which we will use in sections III and V. 



dvii-vy-'-v 
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The sprinkling density S{t) 

Another important quantity for our study is the sprin- 
kling density of events at time t, S{t). An "event" has to 
be understood as a jump from an energy level to another 
one. The sprinkling density of events S{t) is defined as 
the time density distribution of having an event at time t 
-whatever the number of events before t- and given that 
there was one event at t = 0. This is a standard quantity 
defined in the context of renewal theory |8(tI |. For any 
distribution of trapping times ipir) where the trapping 
times are independent random variables -and given that 
it does not depend on the age, i.e on the choice of the 
time origin-, and independently of the type of dynamics 
used, the following formula holds: 



s{t) = m 



dtiS{ti)i>{t - ti) 



(1) 



where ti stands for the time of the last event to have 
taken place before t; the first term i){t) corresponds to 
the special case = 0. 

In the following, we will use the following trapping 
times distribution: iI){t) = j^^^^^pyrfjr , in order to ensure 
that T can take values from to oo. We will also need 
the large time behaviour of the cor resp onding sprinkling 
density S{t). Following the lines of [371, this can be easily 
computed using Laplace transforms. Using the notation 
/(z) to denote the Laplace transform of a function /(i), 
equation (1) is equivalent to 



S{z) = 



^(z) 



1 - V7(z) 

The Laplace transform tpiz) can be computed as: 



dv- 



Two cases have to be considered before taking the limit 
tqz 0. 

If ^ < 1, i){z) = l-r(l-^)(Toz)^ + j^To2 + o(toz). 

Then S{z) ~ j^^J—^^-L^, and for t > tq, S{t) 



with c{fi) — 



sin(7r/i) 



r(i-M)r(M) ^ ■ 
In this regime, S{t) decreases with time; the decrease in 
time of the density of events results in the aging of the 
correlation function Il{tj^,tw + t), which characteristic 
time scale is proportional to the age t^, . 

On the other hand, if > 1, V'(z) = 1+j^toz+o{tqz). 

This imphes that S{z) ~ and that for t > ro, S(t) ~ 

= ^7^. In this non-aging regime, the sprinkling 
density is uniform in time and simply equal to the inverse 
trapping time < r >= dTTip{T) 



Definition of the system of study as a superposition 

of subsystems 

Let us now turn to our system of interest. This 
new system is defined as the superposition of A'' subsys- 
tems, identical to the one introduced previously, each of 
which is defined by the same trapping times distribution: 
"^{t) — . The subsystems are assumed to be in- 

dependent of each other. One can give an interpretation 
of such a model in a real space representation: given a 
system of size L in c? dimensions, we assume that one 
can divide this system into A^ independent subsystems 
of length ^, where A^ = (■|)''. ^ is the typical coherence 
length of the system, and is considered constant during 
the time relaxation. However, we will see that this quan- 
tity is susceptible to evolve during aging. 

During the dynamical evolution, each system relaxes 
independently, and hence contributes to some extent to 
the relaxation of the whole system. We will make the 
following assumptions: (i) all events occuring in a sub- 
system are also defined as individual events for the whole 
system, (ii) all events contribute equally to the relaxation 
of the whole system. 

This translates into the following definitions: 

• (i) The analog for the whole system of iIj{t) will be 
denoted as Pm{t)'. it is the distribution of time intervals 
between all events (i.e trapping times). We will see in 
the next section that this quantity depends in general on 
the age tw\ we will then call it P/v(t, tu;)- 

• (ii) The correlation function of the whole system is 
defined as 



1 ^ 

t) = — Ci {tw , + t) , 



where Ci{tw,tw + ''") is the correlation function of sub- 
system i. 

Before turning to a detailed calculation of P^iTjiw), 
we can invoke an argument of statistics of extremes. If 
T is a trapping time of the whole system, then it seems 
natural to say that r = min{rj}j=i^ where r, is a 
trapping time of each subsystem i. However, this is true 
only if all subsystems undergo one event at some time 
origin, and that one computes the first trapping time of 
the whole system from this time origin. Hence this argu- 
ment definitely excludes aging effects, because it neglects 
any memory effects in the dynamics of the subsystems. 

Having made this approximation, one can follow a 
standard calculation of statistics of extremes, and one 
can find the distribution of time intervals P/v(t): 



PAr(r) = A^V^(t) 



dt'i){t') 



N-l 



(to+t)i+^^ 



To 



By expanding around the most probable value t — 0, 
and setting u = (^Nt/tq, one finds the limiting exponen- 
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e"" for ii <C 1 (i.e t <^ 



tial distribution P(u) 

This corresponds in fact to the convergence of the prob- 
abihty distribution of extremes towards the Weibull dis- 
tribution, in the case where the elementary distribution 
iJj{t) has a finite value for its minimum time t = |38l |. 

In this approximation, the time distribution of events 
of the whole system simply follows a Poisson process, 
with a rate proportional to the number of subsystems. 
In a Poisson process, the conditions of the experiment 
are supposed to remain constant in time, and all events 
are independent of each other. However, in this model, 
although single events in all subsystems are indeed inde- 
pendent of each other, the dynamics is not invariant un- 
der time translations (as can be inferred from the sprin- 
kling density S{t)). As we shall see in the next section, 
this will give rise to more complicated laws for Pjv (''■, tw ) ■ 

Note finally that we have not been able to find a suit- 
able argument of statistics of extremes for C (t^ , + r) 
(the statistics of this quantity in the framework of non- 
equilibrium d yna mics have been related to the Gumbel 
distributions (24l , which are one of the "universal" 
families of probability distributions of extremes). In- 
stead, we will gain information (see section V) by study- 
ing C{t^, tyj + r) as the sum of N random variables, re- 
inforcing the idea that it is not an extremal quantity, but 
rather originates in the contribution of many individual 
events (as was already pointed out in 



III. DISTRIBUTION OF TIME INTERVALS 
BETWEEN ALL EVENTS 

Let Pj^{T,tw) be the probability that an event takes 
place at t„, + T if one took place at tw, in the system 
composed of N independent subsystems. 

In this section, it will be more practical for the com- 
putation to work with the cumulative probability distri- 
bution P^{T,t^) = dr' Pn {t' , tw) ■ By definition, it 
is the probability that the time difference between two 
successive events is larger than r. Similarly, we will use 
Q{t) = /~ dT'tPir') = t^/{tq + tY, which is the prob- 
ability for a trapping time of a subsystem to be larger 
than T, i.e the probability for a subsystem not to change 
trap during the period of time t. 

We now call i the subsystem in which one event has 
taken place at tw Then let {tj}j^i be the {N — 1) times 
of the last events before i„ in the other subsystems j. 
The next event to take place in the whole system will 
either happen in subsystem i or in any other subsystem. 
In order for this next event to occur after a time r, one 
requires the following conditions: 

• (i) subsystem i has to remain trapped between 
and + r, with probability Q{t), and 

'(a) the other subsystems j have to remain trapped 
between tj and tw + t, with probability Q{tw + t) + 
Jjj*" dtjS{tj)Q{tw + T — tj), Q{tio +t) being the contribu- 



tion for the special case tj =0. This last probability is 
in fact equal to Il{tw,tyj + r), the probability for a sys- 
tem to remain trapped between times tw and t^+r (see 
section II). 
Hence, 



pS{T,t^)^Q{T) [n{t^,t^ + T)] 



N-1 



PNir.tw) = 



q{t) [n{t^,t^ + T)] 



N-1 



In the following, we will always consider the case of 
large t^i ^ '''o- We now treat separately two different 
regimes for r: (i) r ~ tq, and (ii) x = j- finite and 
smaller than 1. 



Case T ~ To 



In the regime of interest where ^ to, S{tw) 



1 c(m) 



(see section II). So that the leading term in t„ 



n(i^, t„ + r) ~ 1 - r 5(t„) ~ e-"'5(*"') 
Then, we get the result: 



P^(T,t^) ~ Q(r) e 



(2) 



and 



PN{T,t^) ~ (V.(t) + (iV- l)5(i^)Q(r))e-(^"i)^^('») 

We note that for iV = 1 the result PN{T,t^^) — iP{t) 
is recovered. As N increases, the exponential part in Q 
becomes dominant, introducing a rate that is dependent 
on the waiting time: p = "^^^^^"f!,"^^ . 

This can be interpreted as an effective Poisson process, 
where the number of instances N is replaced (for large 
N) by 



Neff{t^)^N{^ 



In other words, the scaling of l(2Jl suggests that com- 
puting at age t^ the distribution of events of a system 
composed initially of N subsystems is the same as com- 
puting the distribution of events of a "young" system 
composed of Nf,ff{tw) subsystems. In the case of (0, 
one can see that Neff{tw) decreases explicitly with the 
age; in the limiting case where aging disappears (/x — > 1), 
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Nef f is simply a constant equal to N. The idea of a num- 
ber of independent subsystems decreasing with the age in 
non-equilibrium systems is not new. It is intimately re- 
lated to the concept of a growing lengthscale in an aging 
system. If one defines a typical length ^ of a subsystem 
by N — (jY, one has ^eff = LN~^I''', and the depen- 
dence of N^f f on induces the following power-law for 

ee//(i.) = i:iv-V'^(^) ' 

Case X = J- <t:l 

In the regime most accessible experimentally, x — 
^ <C 1, and n(tu,,i^ + t) ~ 1 - c(/.t) ( ^^j , which 
leads to the result: 
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FIG. 1: (Color online) Pf}{T,tm) in the two regimes studied 
for A*' = 1, 2, 5, 10; the values of the parameters are to = 1, 
fj, = 0.5 and = 100. As TV is increased, the departure from 
the power-law is observed. 



and 

(N - 1)c(m)(1 - A*) ^^^^^^^^N-i)cM(^Y-- 

Again the limiting case P^{t, t^) = i^ir) for = 1 is 
recovered. In this regime, when N increases, the dis- 
tribution P^{T,tw) evolves towards a stretched expo- 
nential in T, with a characteristic time proportional to 
the age t^. But as in the case of (jJl, the scaling of 
© suggests that the system at t^^ is equivalent to a 
"young" system composed of Nf.ff{tw) subsystems with 

Let us now be more precise concerning the relevance 
of the quantity Nf,ff{tw)- The physical meaning of the 
effective number of subsystems and of the effective co- 
herence length is the following. Consider the system at 
time ti , with a coherence length ^1 . This system is hence 
made of A^i independent subsystems of size ^1, by defi- 
nition of the model; and each subsystem is a trap model 
that can hop in different states in time, with average hop- 
ping rate S{ti). Each subsystem ages, so that at some 
later time t2, the hopping rate decreases and is equal to 
S{t2)- This also means that per unit time, less subsys- 
tems have hopped than at age ti. In a real space picture, 
the subsystems that hop are then more sparse and far 
away from each other. This induces some kind of en- 
hanced spatial correlation. Then one can make a coarse 
graining of subsystems and define bigger subsystems, of 
size ^2 7 in such a way that the hopping rate of these new 



subsystems is the same as at time ii, which is possible, 
precisely because they are bigger. Hence the total sys- 
tem at time ^2 is now a collection of N2 subsystems of 
size ^2, each of which hops at the same rate as the Ni 
subsystems of size ^1 at time ti. This can actually be 
quantified by a simple argument. The average hopping 
rate of one subsystem at time ti is S{ti). At time t2, 
one defines a coarse-grained subsystem as composed of 
N1/N2 of the former subsystems. Then the average hop- 
ping rate of one coarse-grained subsystem at time t2 is 
S{t2)Ni/N2- The hopping rates are chosen to be equal, 
which leads to A^2 — NiS{t2)/ S{ti). In the special case 
of ti = To and t2 = t^,, one finds: N{t^) = iV(^)i-^, 
which is exactly the relation for Nf^ff(t^) found from the 
previous calculations. 

To conclude this section, the study of the two cases in- 
vestigated above show that for a small system, P^{t, t^) 
(and PN{T,tw)) will still be very close to the power-law 
characterizing one single subsystem. For a very large sys- 
tem, P^(t, t^,) crosses over from an exponential form to 
a stretched exponential form at larger r. 

In general, there will be a crossover in Pn{t, tw) from 
an exponential times a power-law, to a stretched expo- 
nential times a power-law, as r increases. In all cases, the 
distributions become fatter with the age, which allows 
to define an effective number of subsystems Neff{tw) — 

^ (?■) ^ '^^ equivalently an effective coherence length 

For illustration, we plot on Figure 1 the cumulative 
probability distribution P^(t, tu,) in the two regimes 
studied for different values of A^; the values of the pa- 
rameters are tq = 1, = 0.5 and = 100. 
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IV. DISTRIBUTIONS OF TIME INTERVALS 
BETWEEN DECORRELATIONS. APPLICATION 
TO EXPERIMENTALLY ACCESSIBLE DATA 

1. Definition of the quantities of interest 

When one is not able experimentally to identify indi- 
vidual rearrangement events, it may be easier to turn to 
the study of the fluctuations of global quantities such as 
correlation functions. More precisely, in the scattering 
experiments of |l9l |. non averaged correlation functions 
called ci{tw,T) are computed. ci{tw,T) represents the 
degree of correlation between the speckle field scattered 
by the sample at time t^,, and the one scattered at time 

+ T. The time lag t can be given a fixed value during 
the analysis of the data, and one computes the time series 
of c/(tu,,r) as a function of tirne, starting from time 
(see for example Figure 3(a) in [l9j). This allows to com- 
pute the probability distribution of c/, P{ci), for a given 
T and a given age ty^,. This is exactly what we compute 
in section V, if we assume that the function C [t^ , tu, + r) 
of our model can be identified with ci(tw,T). The main 
experimental finding, namely that P{ci) is a negatively 
skewed distribution for a value r small compared to ty^ 
in a non-equilibrium system (a colloidal gel) , is recovered 
in Section V. 

In the analysis of the time series of ci{ty,,T), one can 
also compute the distribution of time intervals between 
significant decorrelations of the system -i.e big downward 
jumps of cj{tw,T)-, which is accessible experimentally. 
Note that these big jumps do not necessarily correspond 
to the same individual "events" studied in Section II and 
section III, and actually correspond to the superposition 
of several of them. However these jumps are the only 
visible manifestation of the individual "events" from an 
experimental point of view, unless new techniques allow 
to visualize in detail and record the spatial rearrange- 
ments of the particles in real space. 

Typically one would like to compute the distribution of 
time intervals between the smallest values of ci{tw, t) -or 
C(ttn, tw + t) in our model-, a threshold value Cth being 



fixed. One can write that C (tw , ^ m i v ~ ^ n 
where n(tu!,tw -\- t) is the number of subsystems that 
changed trap between tw and + t. In this section, we 
will keep the notation A'^ for the number of subsystems, 
since we will always place ourselves at a given age tw ■ The 
infiuence of the dependence of the number of subsystems 
on tw will be discussed in the conclusion and is not crucial 
here. 

The threshold Cth can be chosen such that one se- 
lects only the situations where at least k subsystems have 
changed trap between t^ and t^+T^ so that Cth — ^~ 
and the values of C{tw,tw+''') considered will be less than 
Cth- The distribution of time intervals between succes- 
sive jumps of the correlation of this kind will be denoted 
V'^^^T, tw,T), where T is the time interval variable. 



Before considering the general case, we will first focus 
on two simpler cases: 

• (i) we will first compute T^*-^-* (r, t^;, T), which cor- 
responds to the distribution of times between successive 
decorrelations of maximum intensity (fc = N). We will 
see that the age t^ is not relevant for this quantity, and 
it will be denoted V^^\t,T). 

• (ii) we will also compute V^'''''\T,tw,T), the dis- 
tribution of time intervals between successive decorrela- 
tions of the same intensity, i.e which correspond to the 
case where k subsystems have changed trap between tw 
and tw+T, and k subsystems have changed trap between 
tw + T and tw + T + t, for the first time since tw 



2. Special case of the largest decorrelations 

In this section we want to compute V'-^^{T,tw,T). 
Since we consider only the largest decorrelations, we say 
that an event has occured at tw if C {tw ,tw + t) first 
reaches the value at t„ . We define Vq^^ (r, tw , T) as the 
probability per unit time r that such an event takes place 
at time tw + T, knowing that such an event happened at 
tw This quantity will be helpful in all the following in 
order to calculate V^^^t, tw,T). 

In the case of 'P'^^\t, tw,T), the decorrelations in the 
intervals [tw',tw + t] and [tw + T;tw + T + t] are two 
successive total decorrelations; whereas in the case of 



VQ^\T,tw,T), they may not be successive in time. For- 
mally speaking, 7"^^) (r, tw , T) plays the role of ipiT) and 
r^^\T,tw,T) the role of S(T), where ^p{T) and S{T) 
have been introduced in Section II. 

Throughout the whole section, we will consider the 
regime where t <C T. In this case, a good approximation 
is: 



ri,''\T,tw,T)o,l[TS{T)f 

T 



Therefore, this quantity does not depend on age, and, 
as in section II, one has the relation: 



7'(^)(r,T) =7'(^)(^,T)+ / dt'7'W(T,t')^r^(r,r-t')- 



Equivalently, if z is the Laplace variable conjugated to 
T, the Laplace transforms are related according to 



In the case of interest (/i < 1), we use the fact that 
S{T) ~ c{n)^^ for T > To, so that 
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c(/^)r 



JV 



e 



u.,x—l 



where x = l — N{l—^i), and we have introduced the lower 
cut-off To for the case where the integral is divergent at 
the origin. 

According to the value of N , we will sec below through 
the study of the different cases that V^^\t, T) is a power- 
law with an exponent depending on N. 



(a) Case < x < 1 : N{1 -//)<! 

In this case, the former integral converges when tqz 
0; therefore, 



pW(T,r) ~ 



V \T0 ) 



{T2 + TY+y' 



where T2 = 



Note that the constant in front of the power-law is 
not exact (one would need the expression for all T of 
S{T) to get its exact expression). Moreover, this con- 
stant is smaller than 1, which means that the distribu- 
tion V^^\t,T) is not normalized. This comes from the 
fact that the total number of events (which is equal to 
'Pq^\t,z = 0)) is finite; therefore there is a non-zero 
probability that the time interval between two events is 
infinite, so that V^^^t, T) is not normalized to unity. 



1 



1 



r(l-a;) (Toz)- 



where we have defined Tq = t 



c{h)t 



r(x)r(i-x)- 



Hence, if Tqz < 1, V^^\t,z) ~ 1 - r(l - a;)(To^)^, 
which leads to 



xT§ 



(To+T)i+^' 



Note that in the special case iV = 1, one has x — fj, 
and To = To, and one recovers the power-law with the 
initial exponent 1 + ji. 



(b) Case < j/ < 1 : 1< JV(1 - /u) < 2 

The integral is divergent at the origin in the case where 
X <0, and instead of x, we use for convenience y = —x = 
N{l-^i)-l> 0. 

For < y < 1, or 1 < N{1 - /x) < 2, we use the 
following expansion when tq^ — > 0: 



du 



1 1 



u'^+y y {Toz)y y 



r(i-2/) 



N 



We define the time Ti such that = ^ , and 

we find that P^^^^^) - K^)' " y^^^ ~ 
and 



(c) Case y > 1 : iV(l - > 2 
Now the expansion for tq^; — > of the integral reads: 



f 



du ■ 



1 1 



1 



1 



1 1 



y {Toz)y y~l {Toz)y i yi-y 



r(2-2/). 



in the case where 1 < y < 2. In general, if n = E{y) 
is the integer part of y, the constant part of the expan- 
sion is proportional to T{n + 1 — y) and is followed by a 
(toz)"+^^^ term. Then, following the case 1 < y < 2, 

we find that v!,''\r,z) c i (^)' - ^TfTo^"^. + 

^j^T{2 — y){Tiz)y + o{z'^). In general, there will always 
be a singular term in z*' in between two polynomial terms 
z" and z""*"^. Again in the special case where 1 < y < 2, 



we find: ■P<^^\t,z) 



^(1 - ^^yT. 



2 ~o 



z -h 



^r(2-y)(T2zF + o(z2)). 

As a reminder, let us note that the Laplace transform 
for small 2; of a power-law distribution '<1){t) introduced 
in Section II is, for 1 < < 2: Mz) = 1 rr'^o-z + 

^r(2 - [i){tqzY + o((toz)2). Hence -p(^)(T,z) cannot 
be mapped exactly on this type of function, and one can 
for example add an exponential function to a power-law 
in order to recover the correct expansion up to order z^ . 
This introduces some indeterminacy in the determination 
of V^^\t,T) . However, the behaviour at very large T 
will be: 



7^W(t,z):^ 



i ( nV 



r(i - y) 



y \To J L y \ro ) 

Hence, V^^\t,T) is again a power-law at large T 



pW(T,T)~ 



a 



where a is an indetermined coefficient. This final 
power-law behaviour might be in fact rather hard to ob- 
serve in reality. 
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3. Special case of the decorrelations with same value 



In this section we compute the quantity 
'P^'^''^) {T,tw,T). Like for the previous case, we de- 
fine vl^'^\T,t^,T) corresponding to V'^^'^\T,t^,T), 
which measures the probabihty that k subsystems 
change traps in the interval [t^ + T; + T + r] , knowing 
that k subsystems (not necessarily the same) have 
changed traps in the interval \tw ; + r] . More precisely, 
in this section, we say that an event has occured at 
tw if C{tyj,tw + t) first reaches the value 1 — ;^ at t^,. 
Then Vq^^ (r, tw , T) is the probability per unit time r 
that such an event takes place at time tw + T, knowing 
that an event happened at t^- In this case [k N), all 
quantities will depend on and a relation of the type 
of 1^ is a priori not correct. However, in the limit that 
we consider, where both T and t are small compared to 
tiu, and if k is not too far from N (experimentally, the 
threshold will be taken low enough in order to get rid of 
experimental noise), we will assume that such a relation 
can still be valid. We will write: 



y = n/N, and by replacing J2n=o by N J^'"'' in the limit 
where N is large. 

First, if we consider the lowest order in r only, we 
reduce the expression to: 



V^'''\r,t^,T) = ^ [TSiT)f~- [rSit^+T)Y 

Then by using the Stirling formula A^! ~ 
N^e~^V2TTN, and taking the continuous limit, we 
get: 



(27r)3/2Ari/2 
exp[— A^{(a; — y) ln(x — y) 
+ 2y In y + (1 — a; — J/) ln(l ~ x ~ y) 
- {x-y)\nS{T)-y\nS{T + U,)}] 



dy 



(6) 



By using a saddle point method, we find that the for- 
mer integral is maximized for y* such that 



pf'=)(r,i^,T)=7'('=^'=)(^,i»,T) 



Like in the previous section, we will first com- 
pute v''^''"\T,t^,T) and deduce V^^^^\T,t^,T) from its 
Laplace transform thanks to the relation: 



I \i T I'W, -^J — ~ Ik k) , 



(5) 



In order to compute 'P'^'^\T,tw,T), we call n the 
number of subsystems that changed traps in the in- 
terval + ■''] 7 but stayed trapped in the interval 
[t^ +T,tw + T + T]. Hence k~n subsystems change traps 
both in the intervals [tw ,tw+T] and [tw +T,tw+T+T]; n 
subsystems change traps in the interval [t^+T, t^+T+r] 
but where trapped in the interval [tw,tw + t]; finally, 
N — k — n subsystems were trapped in both the intervals 
[tw ,tw + T] and [tw + T,tw + T + t]. 



Therefore, if k 



sup 



min(/c, N — k), 



{x-y*){l-x-y*)=y*'^{tw,T) 
_ Sir) 



(7) 



where ^{tw,T) = -gfY+iz)- 

We keep the positive solution, having checked that we 
have both y* < x and y* < 1 — x: 



2(7- 



-[-1 + v/l + 4:r(l-a;)(7-l)]. 



Replacing y* in the last expression for 'P'l^'^\T,tw,T) 
leads to a complicated expression that can be simplified 
in our regime of interest, tq <^ T <^ tw 

First, by using the expression for S{T) at large times 
and fi < 1, we find that, for tq <^ T <^ t^, ^{tw,T) ~ 

(^^/^)l-^ and y*{T,tw) ^ (^gj^)'^'- Then, if ^ 
is large enough, we can always assume to have y* 
min(a::, 1 — a;), so that equation © can be expanded in 
y* . Finally, we find the result: 



r^o'-'^ (r, tw,T)^^ [TSiT)f e'^^^^(^) ' ■ 

ZTTT 

The Laplace transform with respect to T reads: 



/^k 

^ Al 



p('='^)(r,t^,T) = ^ ^ [tS{T)] 



k — n 



[1 - r5(T)]" [rSitw + T)]" [1 - r5(t„ + T)] 



N-k-r 



V^'''''\T,tw,z) 



/^k 

27rT 



where 



I (tw ; ^) ; 



Such a finite sum being analytically untractable, we 
will take its continuous limit, by defining x = k/N and 



I{tw,z) 



dT 



X{l-^L)k 



- exp 



K{k) 



T 
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and K{k) = 2^k{N ~k). 

I(tw, z) is a convergent integral for fc(l — /i) < 1, which 
is the first case we consider. 



(a) Case - /i) < 1 

Since we are interested in the case T ^ t^, we will 
make an expansion for zt^^ 3> 1. 

Then we need to consider two separate cases for the 
evaluation oi'p'y^'^) (r, tw , z) from relation (O . We report 
here the main results, and refer to Appendix A for the 
calculations. 



(al) Case < fc < ^it^ 
p('''''^(r, iu,,T) is a power-law: 



AC 1 

r(-/3) T^+p 



(8) 



with a = Mi_£L and /3 = i±i±^ _ i. 



(a2) Case ^ < fc < 

V^^'''\T,t^,T) ~ e^°<iT-,)iP^i*^^^-))^{^)^ (9) 
which is a stretched exponential since cos(y^) < 0. 

(b) Case fc(l - ^l) > 1 

p^^'^) {r^tyj^T) is simply an exponential, with a weak 
dependence on t„ (which vanishes for large t^u): 



4. Generalization 

The most general case, as introduced in the beginning 
of Section IV, consists in computing V'-'^\T,tw,T) the 
distribution of time intervals between successive decorre- 
lations larger than . 

As before, we will use the Vo quantities relative to the 
V quantities we are looking for, and we will assume that 
relations between them similar to equation still hold. 

Therefore, we will start by computing 



E 

k'>k k">k 



(k',k") 



{T,t^,T) (11) 



sub- 
+ t] 



where Vq'' ''^ \T,tw,T) is the probability that k' 
systems have changed traps in the interval [tw , 
and k" subsystems have changed traps in the interval 
[ty^ + T; + T + r]. We refer the reader to Appendix B 
for technical details of the computation, which is similar 
to the one explained in section IV. 3. 
From the calculation, we get 



(''': twj T), 



where k* is a growing function of fc, and is equal to k if 
tS{T) < 1 and fc > N/2. 

Therefore the results of the previous section can be 
used. We conclude that, according to the value of k* 
(which will be equal to k in most cases of interest), 
V^''\T,tyj,T) can either take a power-law form (equa- 
tion((Sl)), a stretched exponential form (equation (jSJ), or 
an exponential form (equation pUI) ). where k has to be 
replaced by fc*. 



r^'^'"'' (r, t^, T) ~ C(i^, r)e-^(*""^)^ (10) 
wtieie L(tuj,T) - ^ —7^ — . 

L U J Tq 

We note here that it is not possible to make a clear 
continuity between the cases k = N and the case fc dif- 
ferent than N treated in this section. This is due to the 
face that we use a saddle-point approximation of I{tw, z), 
which is defined only if K(k) is different from 0, hence 
k is different from TV and 0. In this section, the approx- 
imations are more numerous and are susceptible to be 
valid only for some intermediate values of k. Moreover, 
we are interested in the aging regime and now the correct 
expansion is to be made for zt^ 3> 1, and not anymore 
for zTo <C 1. 



V. PROBABILITY DISTRIBUTION OF THE 
TOTAL CORRELATOR 

In section III, we used the notion of "events" , and cal- 
culated their time distribution, because it was the most 
natural quantity to compute in the framework of the 
model presented here. However, it may be difficult in 
an experiment to identify such events, or even to find a 
reasonable definition of individual events. 

A well-defined quantity that is more readily accessible 
is the distribution probability of a correlator C{ty,,tw + 
r), as has been recently investigated in numerical simu- 
lations [23 |. 

In our model, each subsystem i is a two-level system, 
where the correlation Ci(iu,,tu, -t- r) can only take the 
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values or 1. More precisely, the probability distribution 
of the correlation of a subsystem i is: 



P{C,) = fd{C, - 1) + (1 - f )S{Q), 

where the parameter / coincides with the average value 
of the correlation for the values of t^, and r considered: 
/ Il{tw,tw +t) =< Ci{tw,tw + t) >. In the following 
we will omit temporarily the dependence of / on and 
T, for simplicity in the notations. 

For the whole system (the superposition of the N inde- 
pendent subsystems), the definition of C{tw,tw + t) has 
been given in Section II. One can actually rewrite this 
quantity as: 



m(tv 



N 



N ' 



lu , ito + r) is the number of subsystems that 



where m(f 

remained in the same energy trap between and t^+r, 
and n{twTtw + T) the number of subsystems that changed 
trap between tyj and tyj + t. 

Then, in order to find the value C for the correla- 
tion between times tw and + t, one has to draw m 
subsystems among N instances that are trapped dur- 
ing this time interval, with probability /, and N — m 
subsystems that changed trap during this time interval, 
with probability 1 — /. Hence, the probability distri- 
bution of C, P{C), is simply the binomial distribution: 
P(C) = NP{m), with 



P(m) = C^r(l ~ Z)"^"", 

and Cjy is the binomial coefficient. 

Note that in the limit of a large number of subsys- 
tems TV, the central limit theorem holds and the limiting 
distribution is a Gaussian: 



P{C) 



N 



2^/(1 - /) 



exp 



N{C-ff 
2/(1-/) 



We show on Figure 2 the distributions P{C) for dif- 
ferent values of the average correlation: < C >— 
0.1, 0.5, 0.9. The fuU lines correspond to iV = 10 and the 
dotted lines to iV = 100, showing a closer resemblance to 
a Gaussian distribution for N — 100. 

For comparison with experiments or numerical simu- 
lations, we computed the variance a^{C), the skewness 
s(C) and the kurtosis k{C) of the distribution P{C). 
They can be easily computed by noticing that the cu- 
mulants of independent random variables are additive 
quantities: 




FIG. 2: (Color online) Probability distribution P(C) for dif- 
ferent values of the average correlation: < C >= 0.1 (cir- 
cles), 0.5 (diamonds), 0.9 (stars). The full lines correspond to 
N = 10 and the dotted lines to N = 100. 



1 1-2/ 



J_6/^j-6/-M- 
N /(I-/) 



At this stage, one would like to plot the previous quan- 
tities as a function of t for different values of the age tw , 
in order to compare with existing results in the literature 
^24i22j- For given r and t^j, f = Il{tw,tw + t) is ex- 
plicitly defined using the relation of section II. Moreover, 
we use the results of Section III that suggest that at age 

, the system can actually be considered as a superposi- 
tion of Neff{tyj) — N i^Y'^ independent subsystems. 

Hence in the preceding formulae for (t^(C), s(C), k(C), 
we replace N by N(.ff{tw)- We want to stress that it is 
not proved here that this mapping can be applied to the 
case of the full probability of the correlation functions, it 
is only a phenomenological attempt to try to extend this 
coarse-graining to the calculation of other quantities than 
the distribution of times. This argument actually leads 
to results compatible to what is found in experiments. 

In Figures 3, 4 and 5, cr^(C), s(C), and k(C) are 
shown as a function of r for different values of t^, for 
H = 0.8 and N = 10. As observed in recent experi- 
ments on foams [lil Isij , one can see that the variance 
is maximum at a time r of the order of the relaxation 
time (which in our model is proportional to the age tw), 
and goes to at large t; this maximum increases with 
the age. The skewness is negative at small r, crosses the 
origin at intermediate times when P{C) becomes symet- 
ric and then becomes positive; the negative skewness is 
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more and more pronounced when increases. The kur- 
tosis is positive at small times, then negative, and pos- 
itive again at large r, while the negative part becomes 
more and more pronounced with the age. All these re- 
sults are compatible with recent numerical simulations 
|24| . and seem to be com pat ible with preliminary exper- 
iments on colloidal gels [l9j. Note again that the use 
of Neff(tw) contains crucial information: if N was kept 
constant, the curves would simply be superimposed by a 
T /tw rescaling, in particular, the maxima and minima of 
the variance, skewness and kurtosis would not depend on 
the age. 

Finally, we computed the probability distribution 
P{C) for a given r/t^ ratio - i.e a given value of the 
average correlation < C >-, at different ages, therefore 
for different Neff{tw); the system is made initially of 
N — 1000 independent subsystems. Figure 6 shows such 
probability distributions, for /i = 0.8, / =< C >= 0.8, 
N ^ 1000, plotting P{C)a{C) versus (C- < C >)/cr(C). 
One can see that even for N = 1000, systematic de- 
viations exist, though they may be hard to distinguish 
when an experimental noise is present. Deviations are 
systematic as increases, which is not surprising since 
Neff{tw) decreases with the age; the central limit theo- 
rem is less and less valid as increases. Note that we 
tried to fit these curves with Gumbel distributions like in 
|24l |. but we didn't find a very precise agreement in our 
case. 

Finally, let us note that our results fail to reproduce 
the expected zero limit of the skewness and kurtosis of 
P(C) when —^ oo. These two quantities actually di- 
verge in our case for any finite N, because P{C) tends 
to a delta function S{C) and (j{C) — > 0. However, the 

limit > oo is formally equivalent to Neffit^) — * oo, 

which restores the correct limit because of the central 
limit theorem. This problem would be absent if P(Ci) 
for one individual system had a non-zero variance when 
Ci 0. In practice, this will be always the case, because 
of the presence of white thermal noise. 

VI. DISCUSSION AND CONCLUSION 

We have studied a system made of the superposition 
of N subsystems characterized individually by a slow dy- 
namics and a power-law distribution of times between 
jumps. This system can in turn be characterized by its 
distribution of time intervals between events PAr(r, i^,). 
For small N, the power-law behaviour of a single sub- 
system is the dominant behaviour. For large N, at very 
small time intervals, this distribution is exponential and 
crosses over at larger time intervals to a stretched expo- 
nential. In the two regimes, a natural scaling between N 
and tw appears, which allows to define an age dependent 
effective number of subsystems Neff(tw) which decreases 
with tw, or equivalently, an effective coherence length 



o-o tw-1 




X 



FIG. 3: (Color online) Variance of the probability distribu- 
tion P{C) for fj, = 0.8 and = 10 independent subsystems 
initially, for different values of t^. 
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FIG. 4: (Color online) Skewness of the probability distribu- 
tion P{C) for = 0.8 and A'^ = 10 independent subsystems 
initially, for different values of tw 
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FIG. 5: (Color online) Kurtosis of the probability distribu- 
tion P{C) for n = 0.8 and A'^ = 10 independent subsystems 
initially, for different values of t^. 
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(C-<C>)/a(C) 

FIG. 6: (Color online) Rescaled probability distributions 
P(C) for jj, = 0.8, A'' = 1000 independent subsystems ini- 
tially, < C >— 0.8, and different values of t^n- 

which grows slowly with the age: ^{t^) ^ ^^^'^ . It 
would be interesting to see whether experimental results 
also display the different regimes calculated here. In par- 
ticular, one could observe an exponential or a stretched 
exponential distribution for a young system, and a power- 
law distribution for an aged system. As explained in sec- 
tion III, the interpretation of £,eff{tw) comes from the 
aging of all individual subsystems: the systems that hop 
are more rare when increases, and hence further apart 
from each other; this introduces a relevant lengthscale 
£,{tw) that grows with tw 

We also computed V''''\T,t^,T) the distribution of 
time intervals between successive decorrelations over an 
interval r, larger than This quantity, like P{C), 
is accessible through scattering experiments. Although 
the derivation of this quantity is a bit technical, the re- 
sults can be summarized in the following way. In the 
extreme case where one only counts the largest decor- 
relations (fc = N) , the distribution of time intervals be- 
tween successive decorrelations crosses over from a power 
law with a small exponent at small N (or large tu,) to a 
power-law with a large exponent at large N (or small 
t^). It is hoped that in experiments, these regimes will 
be observed, when one varies the age or other parame- 
ters. For example, if the temperature is decreased in a 
glass, or the density is increased in a gel, or if one in- 
creases the quench rate, one expects that the coherence 
length increases, i.e the number of independent subsys- 
tems decreases. For k N, the results are more difhcult 
to obtain, due to technical complications, and are likely 
to be reliable only for intermediate values of fc. If fc is 
small (i.e one counts only the decorrelations correspond- 
ing to a few subsystems), the distribution V^'''' {T,t^,T) 
is a power-law; if k is larger, it is a stretched exponential 
and for k large it is an exponential. Equivalently, since 
one expects the number of subsystems to decrease with 



the age, 'P*^'^-' (t, T) is typically an exponential at small 
tu,, then a stretched exponential and finally a power-law 
at large t^^ , for a given value of the correlation threshold 
Cth^l- k/N. 

Finally, we also computed analytically the probabil- 
ity distribution of the correlator P{C), which reproduces 
many features of the experimental findings of |0| and 
of the numerical results of We used the fact (not 

shown rigorously) that at tw the system can be de- 
scribed as a collection of Neff{tw) independent subsys- 
tems. No rescaling to some time-evolving Gumbell distri- 
bution seems to be relevant here. Since the full probabil- 
ity is calculated here, it would be interesting to compare 
it more quantitatively to experimental results. In partic- 
ular, we find that the variance grows with the age like 

^ ^^eff ■ 

Although these results are derived from a simplified 
theoretical model, we think that such a description cap- 
tures the essential intermittent character of the relax- 
ation of a realistic glassy system, together with the ex- 
istence of a coherence lengthscale growing with the age, 
due to the decrease of the rate of events in time. We 
hope that these findings will be helpful in analyzing more 
quantitatively some experimental results. In particular, 
the results of Section IV and Section V might be of inter- 
est for the Time Resolved Correlation experiments of . 
However, since the individual trap model is probably a 
very crude model for the relaxation of realistic materials 
(< C(tyi,,tw + r > is a priori unlikely to coincide with 
the experimental average relaxation function of colloidal 
gels), it may be interesting in a further study to replace 
it with a more realistic dynamical model for jammed col- 
loidal gels 40|. Concerning the intermittency dynamics 
studied in 20], the quantities calculated in Section III 
are maybe the most relevant ones to compare with the 
experimental results. 

The model studied here has the major disadvantage 
that the relevant coherence length of the system is intro- 
duced by hand through the superposition of several sub- 
systems. However, we were able to extract useful infor- 
mation from this picture. It would be more satisfactory 
to solve a model where an aging coherence length builds 
up during the dynamics, resulting from microscopic dy- 
namics. Candidates could be special cases of kinetically 
constrained models, though analytical calculations will 
be quite hard [4]| . However, numerical simulations pro- 
vide now a useful alternative tool for the precise investi- 
gation, in microscopic models of glasses, of the quantities 
calculated in this paper ,24] . 
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APPENDIX A 

(a) Case fc(l - fi) < 1 

In this case, we can perform a saddle-point calculation 
in order to evaluate z), which leads to: 



with A 

— l-M 



[(1 - f,)K{k)/2]-'^, B = [{1 - 
- MWO!, p ^ _ 1, and 



Since we are interested in the case T <C , we expand 



a J) 



1 



B 



I{tyj,z) for ztw > 1: Iitw,z) 
Then we need to consider two separate cases for the eval- 
uation of 'p(''^*=)(r, t^, z) from relation Q. 



(al) Case < fc < ^^tii. 



(b) Case fc(l - > 1 

Finally, we consider the case fc(l — /i) > 1: then we 
have to introduce the lower-cut-off tq in which 
leads to: 



du 



The saddle point solution is u* — [(1 — 
/i)i^(fc)/2]T+f' (zt^)TT?r. In the limit of very large 
ztw, one will at some point reach the situation 
where u* < tqz. Hence, for zt^^ ^ 1, the integral 
will be best evaluated by its lower cut-off value: 

. Since T ^ tq, 



I{t^,z) i^e-o^e^«(^^)" 



we 



finally get: 'P^'''''\T,t^, z) ~ c(t^^ 



re 



2irT 



Finally, P^'^''^^ (r, t^, T) is simply an exponential, with 
a weak dependence on (which vanishes for large tw)'- 



APPENDIX B 

For the computation of 'p'^ \T,t^,T), we will 
proceed in the same way as in the last section for 

Ik k) f 

Vq ' (T,t^,T). Using again combinatorial arguments, 
we can write that: 



For < fc < 2^13^ (i-'^ P < 0), one gets for zi^ > 1: 



c(p)t 



Therefore, V''^'^^{t, t,,,T) is a power-law: 



fc' 



[1 - tS{T)]" [TS{t^ + [1 „ rS{t^ + T)f-''"- 



ZTTT 



Af?,, 



r(-/3) T^+f^ 



(a2) Case ^ < fc < 



where Uinf — sup(0, fc' — fc") and Ugup = inf (fc', N — fc"). 

As before, we keep only the leading order in r and go 
to the large N limit by introducing x' = x" ~ and 
y = j^. This leads to: 



In this case, we find that V^^''^'' {T,t^, z) 

ik 



exp 



K.(t„,T) 



, with /C(t„,r) = 



The evaluation of the inverse Laplace transform leads 



to: 



^fc"-i 



(27r)3/2iVi/2 



dy e 



-NS[x\x" ^y) 



where 



/3 



which is a stretched exponential since cos(y^) < 0. 



S{x' ,x" ,y) = y\ny+ {x" - x' + y) \n{x" - x' + y) 
+ {x' - y) ln{x' ~y) + {l~ x" - y) ln(l - x" - y) 
- {x' - y) In S{T) - {x" ~ x' + y) In S{t^, + T) . 



14 



We use a saddle-point approximation to minimize 
Six',x",y): If Ij,. =0 leads to: 



{x' - y*){l - x" ~ y*) = y*{x" - x' + y)^{t^,T) 

which is the generalization of equation iQ. 
The positive solution is therefore: 



^ 2(7-1 
+ v/[l + (7-l)( 



l-(7-l)(x"-a;') 



+ 4a;'(l - a;")(7 - 1)] 



In 



1 -a;'^ 



1 - 2x" 



Vx'*(l-.T'*)7(t„,T) 



ln[r5(T)]. 



If tS{T) = 1, this equation has the only solution 



i. For tS{T) > 1, there is a solution 



a;'* > i 



and for tS{T) < 1, there is a solution x'* < ^. In gen- 
eral, X* — sup(a;,a::'*) will maximize the expression for 

'Pq^\t, tyj,T), and can be estimated numerically. In par- 
ticular, if tS{T) < 1 and k > N/2, we have x* = x. 
Finally, 



and it is easy to check that Xinf < y* < Xsup ■ Note that 
in this equation and in the following, we drop the (t^, T) 
dependence of "/(tj^,T) for convenience. 

We now come back to the calculation of Vq'''^ (r, ty, , T) 
through equation Hll|) . For large N, we have: 



dx' 



where k* = sup(/c, k' 



dx"V^i'-''"\T,U 



Direct integration is not possible here, so we use again 
a saddle point treatment of the double integral. The sad- 
dle points x'* and x"* are solutions of the set of equations 
formed by: §§t\{x",x"') = and ^\(x",x"') = 0, or: 



or 

dx' 



In 



y*{x" 



{1 - x" - y*){x' - y*) 



In 



y*) 



{x' - y*) 



(12) 



In 



TSiT + t^){l-x" 



9y* -^^ y*{x" -x' + y*)^ 
dx" ^ (1 - x" - y*){x' - y*) 

"(13) 

Here, ^p-\(x'* ,x"') and ^ijt\(x" ,x"') can be calculated 
using the expression for y* . 

At this point, taking the ^ 1 limit enables us to 
make simplifications in equations (|12|l and H13|l . A rather 
lengthy calculation leads to the conclusion that there is 
no couple of solutions {x'* , x"* ) compatible with equa- 
tions l(T^ and ((T^ . and such that x'* ^ x"* . 

Consequently, we keep only the terms such that k' = 
k" in the sum of (|ll|l . and try to find the value for 
k' that maximizes this sum, knowing the expression of 
v''Q'''''\T,t^,T) from section IV.3: 



V^''\T,t^,T) (X I da;' exp{iV[-x'lna;' 



- (1 - x') ln(l - x') + x' \nS[T) + 2(^^iL^) V^]} 
The saddle-point equation for x'* is now: 
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